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In 1959 Weibel demonstrated that when a QED plasma has a temperature 
anisotropy there exist unstable transverse magnetic excitations which grow ex- 
ponentially fast. In this paper we will review how to determine the growth rates 
for these unstable modes in the weak-coupling and ultrarelativistic limits in which 
the collective behavior is describable in terms are so-called "hard-loops". We will 
show that in this limit QCD is subject to instabilities which are analogous to the 
Weibel instability in QED. The presence of such instabilities dominates the early 
time evolution of a highly anisotropic plasma; however, at longer times it is ex- 
pected that these instabilities will saturate (condense). We will discuss how the 
presence of non-linear interactions between the gluons complicates the determina- 
tion of the saturated state. In order to discuss this we present the generalization 
of the Braaten-Pisarski isotropic hard-thermal-loop effective action to a system 
with a temperature anisotropy in the parton distribution functions. The resulting 
hard-loop effective action can be used to determine the time and energy scales as- 
sociated with the possible saturation (condensation) of the gluonic modes. We will 
also discuss the effects of anisotropics on observables, in particular on the heavy 
quark energy loss. 



1. Introduction 

In the last few years a more or less standard picture of the early stages of 
a relativistic heavy ion collision has emerged. In its most simplified form 
there are three assumptions: (1) that the system is boost invariant along the 
beam direction, (2) that it is homogenous in the directions perpendicular to 
the beam direction, and (3) that the physics at early times is dominated by 
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gluons with momentum at a "saturation" scale Qs which have occupation 
numbers of order l/a^. The first two assumptions are reasonable for de- 
scribing the central rapidity region in relativistic heavy-ion collisions. The 
third assumption relies on the presence of gluonic "saturation" of the nu- 
clear wavcfunction at very small values of the Bjorkcn variable a; In this 
regime one can determine the growth of the gluon distribution by requir- 
ing that the cross section for deep inelastic scattering at fixed does not 
violate unitarity bounds. As a result the gluon distribution function satu- 
rates at a scale Qs changing from 1/fc^ \og{Ql/k'\)/as. Luckily, despite 
this saturation, due to the factor oil/ as in the second scaling relation the 
occupation number of small-a; gluonic modes in the nuclear wavefunction 
is still large enough to determine their distribution function analytically 
using classical nonlinear field theory ^. 

These distribution functions are used as input for the subsequent ther- 
malization of the quark-gluon plasma. At short times, tq ^ Q^^i the gluon 
distribution function is isotropic, however, after the system begins to ex- 
pand it rapidly develops a large anisotropy between the transverse and 
beam directions with p±_ due to the fact that initially ~ r~^. In 

the weak-coupling limit the assumptions above have been used by Baier 
et al. in an attempt to systematically describe the early stages of quark- 
gluon plasma evolution in a framework called "bottom-up" thermalization 
^. Using coUisional mechanisms they are able to show that during the first 
stage of evolution hard gluons scatter out-of-plane counteracting the effect 
of the expansion of the system reducing the rate at which the longitudinal 
momentum decreases instead to Pz ~ qV^t~^/^. Although less extreme 
in terms of the rate at which the longitudinal momenta decreases this sce- 
nario still results in considerable momentum-space anisotropics in the gluon 
distribution function. In such anisotropic systems it has been shown that 
the physics of the QCD collective modes changes dramatically compared 
to the isotropic case and instabilities are present which can accelerate the 
thermalization and isotropization of the plasma 3,4,5,6,8,9,io,ii_ 

In addition to their role in isotropization and thermalization it is also 
important to have an understanding of what the impact of such instabili- 
ties will be on observablcs at RHIC and LHC energies. Here wc will discuss 
one test observable - the heavy fermion energy loss 1^,13^ Qj-^g niight be 
worried that there is a fundamental problem with perturbation theory since 
the presence of instabilities naively causes the calculation of the soft con- 
tribution to the heavy quark energy loss to be divergent; however, there 
is protection mechanism dubbed "dynamical shielding" which renders the 
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collisional energy loss finite for QED and QCD i2,i3 However, at realistic 
couplings the presence of instabilities and associated poles on neighboring 
Riemann sheets ^'^ causes significant changes in the soft energy loss contri- 
bution for both QED and QCD. In fact, as we will discuss, the unphysical 
polos can cvcm change the sign of the heavy quark energy loss at low mo- 
mentum turning it instead into energy gain. 

2. Anisotropic Gluon Polarization Tensor 

We consider a quark-gluon plasma with a parton distribution function 
which is decomposed as 

/(p) = 2Nf (n(p) + n(p)) + 4Ncng{p) , (1) 

where n, n, and Ug are the distribution functions of quarks, anti-quarks, and 
gluons, respectively, and the numerical coeSicients collect all appropriate 
symmetry factors. Using the result of Ref. [6] the spacelike components 
of the high-temperature gluon self-energy for gluons with soft momentum 
(fc ~ gT) can be written as 

where the parton distribution function /(p) is, in principle, completely 
arbitrary. In what follows we will assume that /(p) can be obtained from 
an isotropic distribution function by the rescaling of only one direction in 
momentum space. In practice this means that, given any isotropic parton 
distribution function fiso{p), we can construct an anisotropic version by 
changing the argument of the isotropic distribution function 

/(p) = yr+e/iso (vp'+^(p-fi)2) , (3) 

where the factor of vT+T is a normalization constant which ensures that 
the same parton density is achieved regardless of the anisotropy introduced, 
n is the direction of the anisotropy, and ^ > — 1 is an adjustable anisotropy 
parameter with ^ = corresponding to the isotropic case. Here we will 
concentrate on ^ > which corresponds to a contraction of the distribution 
along the fi direction since this is the configuration relevant for heavy-ion 
collisions at early times, namely two hot transverse directions and one cold 
longitudinal direction. 
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Making a change of variables in (2) it is possible to integrate out the 
|p I -dependence giving ^ 

where cos6'„ = k • n and /j,^ = vT+T "^i) > 0- The isotropic Debye mass, 
mo, depends on /iso- In the case of pure-gauge QCD with an equilibrium 
/iso we have = gT. 

The next task is to construct a tensor basis for the spacelike components 
of the gluon self-energy and propagator. We therefore need a basis for 
symmetric 3-tensors which depend on a fixed anisotropy 3-vector n' with 

= 1. This can be achieved with the following four component tensor 
basis: A'^ = S^^-k'k^/P, B'^ = k^k'/k^, C = fi'n^ jh^, and D^^ = k'W + 
k^h^ with = A^^h^. Using this basis we can decompose the self-energy 
into four structure functions a, (3, 7, and 5asn = aA-|-/3B-|-7C-|-(5D. 



3. Collective Modes 

As shown in Ref. [6] this tensor basis allows us to express the propagator 
in terms of the following three functions 

A-'{K)=k''-u;^ + a, 

where 2^1 = ± VTF^^4((^rr7TF);5^rp^2j2y and = a + /3 + 
7 + k^. 

Taking the static limit of these three propagators we find that there arc 
three mass scales: m± and TOq. In the isotropic limit, ^ ^ 0, = m?_ = 
and m\ = m%. However, for ^ > we find that < for all ^ Tr/2 
and rn?_ < for all |0„| < 7r/4. Note also that for ^ > both and m?_ 
have there largest negative values at 6'„ = where they are equal. 

The fact that for ^ > both and can be negative indicates 
that the system is unstable to both magnetic and electric fluctuations with 
the fastest growing modes fociiscd along the bcamlinc (f?„ = 0). In fact it 
can be shown that there are two purely imaginary solutions to each of the 
dispersions relations A^^{K) = and AZ^{K) = with the solutions in 
the upper half plane corresponding to unstable modes. We can determine 



''Explicit analytic integral expressions for a, fS, 7, and S can be found in Ref. [6]. Addi- 
tionally, analytic expressions in the small-^ limit and for propagation along the anisotropy 
direction can be found in Refs. [6] and [10], respectively. 
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the growth rate for these unstable modes by taking w 
the resulting dispersion relations for r(A;). 

In Fig. 1 we plot the in- 
stability growth rates, and 
F_, as a function of wave 
number for ^ = 10 and 6'„ = 
n/8. Note that both growth 
rates vanish at fc = and 
have a maximum F* ^ A*/10 
at fc* ~ /z/3. The fact that 
they have a maximum means 
that at early times the sys- 
tem will be dominated by un- 
stable modes with spatial fre- 
quency l/fc*. 

4. The Hard-Loop Effective Action 



iT and then solving 



0.04 




Fig. 1: Instability growth rates. 



In the previous section by only considering the polarization tensor we im- 
plicitly used a linear analysis which suffices only during very early times. 

As the instabilities grow non-linear interactions become more important 
and can halt their growth. In QED simulations such a saturated state 
does emerge and is called a Bernstein-Greene-Kruskal wave; however, in 
QCD any corresponding saturated state would necessarily be a much more 
complex beast. Answering these q\.icstions requires knowledge of the full 
hard-loop effective Lagrangian in an anisotropic system. Luckily, it is possi- 
ble to obtain a simple expression which generates all hard- loop (HL) vertex 
functions ^ 



1 

. - F F^" 4- — 
4 ^ 2 



/(P) 



,p iPi ^ • {p-Dy 

For example, from this we can obtain the hard-loop 3-gluon vertex 



(5) 



(6) 



/p ^P-qp-r p-kp-q^ 

In contrast to the isotropic n-gluon HL vertices which vanish in the static 
limit, the anisotropic n-gluon HL vertices can be non-zero. This opens up 
the possibility for additional saturation mechanisms which are governed by 
soft non-abelian physics sensitive to these higher vertices. 

Note that Arnold and Lenaghan have shown that if one considers fluc- 
tuations in the vector potential which are directed along the anisotropy 
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direction and only depend on the longitudinal coordinate then it is possible 
to reduce the effective lagrangian (5) to a quadratic form which simplifies 
the analysis considerably They then make some additional assumptions 
including ignoring time-depedence in the self-energy in order to construct a 
toy-model which is then used to study the subsequent evolution of instabili- 
ties in the system. Their results indicate that in non-abelian gauge theories 
that possible non-abelian saturated states are metastable and as a result 
instability growth forces the system along abelian directions. However, it 
should be noted that by ignoring the time dependence in the potential 
as they have done, the instability growth rate in their model is finite at 
k = 0. This is not the case when time dependence is included in the full 
analysis and instead one finds that the growth rate at fc = then vanishes 
as shown in Fig. 1. This difference could have a significant impact on the 
abelianization particularly on the "abelianization length" measured in their 
simulations. 

5. Effects on observables 

Since the unstable modes manifest themselves as poles of the propagator in 
the static limit, it has been argued that the presence of these instabilities 
in general prohibits the calculation of observables in a perturbative frame- 
work, since those quantities would be plagued by unregulated divergences. 

However, it turns out that at least one observable, namely the collisional 
energy loss of a heavy fermion, is protected from these unregulated diver- 
gencies by a mechanism dubbed "dynamical shielding" (see Ref. [7] for 
a proof in the case of ^ ^ 1 and ^ ^ oo). As a consequence of dynamical 
shielding (which incidently is somewhat similar to dynamical screening of 
the magnetic interaction in QCD, hence the name), the collisional energy 
loss becomes calculable pcrturbativcly also in the anisotropic case. There- 
fore, we will use the collisional energy loss as "test observable" to learn 
something about the effects of the instabilities on physical quantities. 

Schematically, the soft contribution to the collisional energy loss is given 

by 

where A{z,k,6n) = v'^A^^v^ is the relevant contraction of the propagator 
and z = li-v with v the velocity of the heavy fermion. The instabilities 
would in principle affect the integrand for Re z = causing singularties to 
be encountered along the momentum integration path. In order to see that 
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Fig. 2: Sketch of the complex z plane including the extension of the logarithm 
to the unphysical sheet. Also shown is how a pole in the unphysical region 
(mountain) has effects felt on the physical sheet. The black line indicates 
where the two sheets are joined together. 

there are, in fact, no singularities we focus on terms of the form (g^ + 
Due to the fact that in the static hmit the structure function a is negative- 
valued we could encounter a singulartity at = a; however, if one takes 
the static limit of a carefully we find 

\\-ma{z)^M'^(^-\^iDz) + 0{z^) , (8) 

z— »0 

where M and D depend on the angle of propagation with respect to the 
anisotropy vector and the strength of the anisotropy. As long as D is non- 
vanishing the singularities are regularized because of the z in the numerator 
of Eq. (7) and we call the singularity "dynamically shielded" . This effect 
can be shown to regulate all singularities which would have naively come 
from the presence of instabilities. Therefore, there is no catastrophic effect 
from the instabilities on the collisional energy loss. However, there is an 
effect which is associated with the presence of the instabilities; in order to 
understand this effect it is however necessary to extend the analysis of the 
collective modes of an anisotropic plasma to arbitrary Riemann sheets. 

The extension to the unphysical Riemann sheets is achieved by continu- 
ing analytically the structure functions beyond their cut structures (which 
for finite f is a logarithmic cut running along the real z axis for < 1, in 
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complete analogy to the isotropic case, but a square-root cut for < sin^ On 
for ^ oo)^". More precisely, since one can continue the structure functions 
from either below or above the cut, the physical sheet has two neighboring 
unphysical Riemann sheets, which - at least in the case of finite ^ - could 
themscilvcs be extended to reproduce the well-known "spiral staircase" form 
of the complex logarithm. Once the structure functions on the unphysical 
sheets are known, one may conduct the analysis of the collective modes as 
was done on the physical sheet, finding - among others - singularities for 
complex z on the neighboring unphysical sheets that extend to the region 
of spacelike momenta 

But why bother about these modes, given that they do not "live" on the 
physical sheet? To see that there can be an effect on physical quantities, 
imagine once again the spiral staircase spanned by the complex logarithm: 
the physical sheet would correspond to the region covered by the spiral 
plane from the ground floor to the flrst floor, while the unphysical sheet 
where the extra quasiparticle mode lives would correspond to the region 
first to second floor. Since the existence of such a mode corresponds to a 
pole in the propagator, there is has a mountainous dent (singularity) in the 
spiral plane somewhere from the first to the second floor. However, since 
the mountain has a finite width, its base can be felt also below the second 
fioor, especially if its peak is near the first floor (see Fig. 2 for a sketch). 
Therefore, the nearer a mode on a neighboring unphysical sheet comes to 
the border of the physical sheet, the more will it affect the propagator on 
the physical sheet and consequently any physical observable sensitive to the 
relevant region of phase-space (in this case spacelike momentum). 

This is precisely the situation encountered for the collisional energy loss: 
there are unstable modes on the physical sheet which, once the momentum 
k becomes larger than some critical value (see Fig. 1), move onto the neigh- 
boring unphysical sheets and subsequently influence the propagator on the 
physical sheet by the mechanism explained above. In fact, it turns out that 
the the contributions from these unphysical modes may drive the leading- 
order perturbative results for the collisional energy loss negative for small 
fermion velocities, corresponding to an energy transfer from hard to soft 
scales. Note, however, that this negative energy loss is not due to the heavy 
fermion in question having a sub-thermal velocity but is instead connected 
to the instabilities of an anisotropic system, as already noted in Ref. [15]. 

For further discussion of the heavy quark energy loss in an anisotropic 
quark-gluon plasma we refer the reader to Ref. [13]. In addition to eluci- 
dating the impact of unstable modes and associated poles on unphysical 
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Riemann sheets, in that paper we showed that for anisotropic systems the 
heavy quark colhsional energy loss has an angular dependence which in- 
creases as the coupling and/or anisotropy is increased. Quantitatively, for 
as = 0.3 and a 20 GeV bottom quark we found that the deviations from 
the isotropic result were on the order of 10% for ^ = 1 and of the order of 
20% for ^ > 10. When translated into the difference between longitudinal 
and transverse energy loss this results in a 10% difference at ^ = 1, a 30% 
difference at ^ = 10, and a 50% difference at ^ = oo. 

6. Conclusions 

In this paper we have attempted to summarize the developments which 
have occured in the last two years regarding the physics of a quark-gluon 
plasma which has a momentum-space anisotropy in the underlying parton 
distribution functions. We have presented expressions for the self-energy 
and full effective action in this case and discussed the emergence of un- 
stable modes which cause exponential growth of the gluon field along the 
anisotropy direction. In addition, we discussed the expected effects such 
anisotropics have on observables concentrating on the heavy quark energy 
loss test observable. 

As we have tried to illustrate here the physics of anisotropic plasmas is 
qualitatively different from isotropic ones. In the weak-coupling limit the 
short-time behaviour of anisotropic plasmas is dominated by the; growth 
of unstable modes. In order to properly understand thermalization and 
isotropization of a quark-gluon plasma it seems necessary to take into ac- 
count the effect of such unstable modes from the earliest times. We note in 
closing if a Vlasov description is applicable at very early times (tq ^ Qj^) 
then instability broadening of the longitudinal momentum is much more 
effective than collisional broadening and it is possible that anisotropics 
due to expansion of the system are, in fact, never generated or reduced 
significantly. However, a more conservative approach calls for the Vlasov 
description to only be applied once the hard gluon occupation number be- 
comes less than one which according to collisional arguments occurs at 
a time ri ^ a~^^'^Qj^ . If this more conservative approach is applicable 
then the anisotropy at ti can be parametrically estimated to be ^ ~ 

^ a7^^^ if collisional broadening is taken in account) which results in 
a strongly anisotropic system in the weak-coupling limit. Whichever ap- 
proach is used, however, it is clear that it is necessary to take into account 
the unstable modes. 
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